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Employing the interband exciton scattering model, we have performed a numerical study of the 
direct photogeneration and population relaxation processes contributing to carrier multiplication 
(CM) in nanocrystalline and bulk PbSe. We argue that in both cases the impact ionization is the 
main mechanisms of CM. This explains weak contribution of the direct photogeneration to the total 
quantum efficiency (QE). Investigated size scaling of QE in NCs and comparison to the bulk limit 
provide microscopic insight in the experimentally observed trends. 



Recent demand in new efficient photovoltaic devices 
has turned significant attention to the problem of car- 
rier multiplication (CM) in semiconductor materials [TJ 
[2]. Early reported ultrafast spectroscopic studies of 
nanocrystals (NCs) [3H5] suggested that nanostructur- 
ing is a direct way to highly increase the CM quantum 
efficiency (QE), i.e., number of excitons produced per 
absorbed photon. This was supported by the following 
arguments: An enhancement of the Coulomb interaction 
between carriers due to their spatial confinement should 
lead to more efficient biexciton production. Relaxation 
of the quasimomentum conservation constraint by break- 
ing the translational symmetry should open additional 
pathways for CM while the intraband cooling is slowed 
down by the presence of the phonon bottleneck. More 
recent reports claim much lower QE or the absence of 
CM in NCs [6H9] . The controversy possibly rise from ex- 
perimental inaccuracy |5], sample-to-sample variation in 
surface properties [10] , and contribution from extraneous 
effects, e.g., photocharging [TTJ[T2]. Comparison to the 
bulk PbS and PbSe spectroscopic measurements shows 
that bulk QE exceeds validated QE in NCs if compared 
on the scale of absolute energy [13J. These issues call for 
the reassessment of the quantum confinement role in the 
CM dynamics [7]. 

From theory point of view, multiexciton photogener- 
ation [5j [14] , coherent multiexciton production [TSJ [16] , 
and incoherent impact ionization (II) [T7lfT9] , are cur- 
rently under debate as primary mechanisms of CM. The 
ab initio calculations on small clusters (< 1 nm) demon- 
strated the role of strong Coulomb correlations in multi- 
exciton photogeneration and emphasized the role of fast 
exciton-biexciton dephasing [20]. Atomistic calculations 
focused on small diameter (< 3 nm) NCs show that CM 
is dominated by II processes [I9j l2T] , An extrapola- 
tion of these calculations to larger NCs and comparison 
to the bulk have been reported recently [22]. Effective 
mass models have been used to study direct photogener- 
ation [23. and coherent CM dynamics in larger (4 — 6 nm) 
diameter NCs. Reported experimental studies consider 
even larger diameter (< 10 nm) NCs. However, no sys- 
tematic investigation of all predicted CM pathways con- 



tributing to the ultrafast spectroscopic signal (i.e., pho- 
togeneration and population relaxation) in various size 
NCs and comparison to the bulk limit has been reported 
yet. 

In this letter, we report on such study employing the 
interband exciton scattering model (IESM) [24] param- 
eterized for PbSe materials. The IESM treats on equal 
footing the biexciton photogeneration during interaction 
with the pump pulse and subsequent population relax- 
ation processes. This allows us to evaluate the contri- 
butions of all predicted pathways to the total QE. To 
clarify the quantum confinement effect on the key quan- 
tities determining QE, comparison with the bulk limit is 
reported on the absolute energy scale [7JE2]- In general, 
performance of the photovoltaic devices can be charac- 
terized by variously defined power efficiency [I2j [22] , in 
particular using the photon energy scale normalized per 
NC/bulk band gap energy [25 . Since we focus on funda- 
mental mechanisms of CM rather than applications, the 
latter unitless scale is not used below. 

We consider an ensemble of NCs in which no more than 
one photon can be absorbed per NC leading to produc- 
tion of either exciton or biexciton state. Hence, QE can 
be defined as 

QE = (2N XX + N X )/(N XX + 7\y , (1) 

where N x and N xx are the ensemble non-equilibrium ex- 
citon and biexciton populations, respectively. For an en- 
semble of NCs, the bulk limit can be defined as the ther- 
modynamic limit: V — » oo, V/v — » oo, and v = const. 
Here, V is a NC volume, v is the unit cell volume, and 
ratio V/v gives the number of unit cells. Since, QE 
(Eq. 0) is volume independent in the bulk limit, it is re- 
ferred to as the intensive parameter. As a result, its scal- 
ing with NC size compared to the bulk limiting value pro- 
vides convenient measure of the quantum size effects. In 
general, the populations N x and N xx as well as DOS and 
Coulomb interactions (determining N x and N xx ) are not 
intensive quantities. To make proper comparison with 
the bulk limit, below, we eliminate their volume scaling 
by multiplying them with powers of (v/V). 
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We start by defining the populations n x and n 
of a-th exciton and k-th biexciton states, respectively. 
Their time evolution during the interaction with the 
pump pulse (photogeneration) and subsequent popu- 
lation relaxation are numerically calculated using the 
weak Coulomb limit of IESM [24 . To account for high 
DOS region in which CM dynamics occur, we recast n x , 
and n k x to the quasicontinuous frequency representation: 
n x{u) = E a >i K5(u-uZ) and n xx (u) = E fc >i n x k x 5(u- 
oo k x ). Accordingly, the intensive populations entering 
Eq. ([!]) become defined as N x = (v/V) J °° duj n x (uo) and 
N xx = (v/V) J °° du n xx (u). 

In most ultrafast CM experiments, the pump duration 
exceeds typical dephasing time leading to the continuous 
wave (CW) excitation regime. Our simulations show that 
in agreement with Ref. [16], high DOS leads to the sup- 
pression of the oscillating terms that contribute to n x and 
n xx _ rpkg iesm predicts interference of the photogener- 
ation pathways contributing to n k x [24]. However, the 
simulations show that due to high DOS, the sign- varying 
terms cancel out fully eliminating the interference. Tak- 
ing all this into account, the CW photogenerated biexci- 
ton population can be represented as 
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and the leading contribution to the exciton population 
as N x (uj prn ) = Ap x (ujpm)- Here, A is proportional 
to the pump fluence, u prn is the pump frequency, and 
7 is the dephasing rate between exciton and biexciton 
states. p x (uj) = (v/V) 2 ^ a 5(u - u x ) and p xx (u) = 
(v/V) 4 S(uj — uJ k x ) are intensive exciton and biexciton 
DOS, respectively. The associated optically allowed ex- 
citon DOS is p x (uj) = (v/V) J2 a \p x \ 2 5(uj — uj x ) and joint 
biexciton DOS is p xx (uj u uj 2 ) = (v/V) 4 ^2ki \Pki\ 2s ( u i 
uu k x )5(uu2 —uo xx ). They depend on the interband exciton, 
p x , and the intraband biexciton, p k f, transition dipoles. 

Key intensive quantity in Eq. ([2| is the effective 
Coulomb interaction which we define as the r.m.s. of 
the interband Coulomb matrix elements connecting the 
states with frequencies falling into [uji,uji + duj\) and 
[o;2,^2 + du)2] intervals, 
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quency, and V x £ x is the Coulomb matrix element be- 
tween these states. The size scaling of V x ^ x in NCs 
reflects the net result of the Coulomb matrix elements 
scaling, relaxation of the momentum conservation con- 
straints and the appearance of new selection rules associ- 
ated with the confinment potential symmetry. In Eq. ([2| , 
V x j?f(uj) denotes effective Coulomb coupling between the 
ground and biexciton states. 

Three terms in the r.h.s. of Eq. ^ describe the associ- 
ated photogeneration pathways [24] : In the first term the 
optically allowed exciton DOS, p x (oo vrn ), depends on the 
pump frequency indicating resonant production of single 
excitons subsequently scattered to the biexciton states. 
This is the indirect biexciton photogeneration pathway 
predicted by IESM [24]. The second term contains biex- 
citon DOS, p xx (oj p m), depending on the pump frequency 
indicating the direct biexciton photogeneration through 
virtual exciton states. Finally, the last term describes 
the pathway involving direct biexciton generation by op- 
tical stabilization of the scattering processes between the 
ground state and biexciton manifold. The last two path- 
ways have been first considered in Refs. [5] and [14] . re- 
spectively. 

The population relaxation dynamics is described by a 
set of kinetic equations introduced in Ref. [24 and fur- 
ther transformed to the quasicontinuous frequency rep- 
resentation. They contain II and Auger recombination 
rates 
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Here, uo x (uo^ ) is a-th exciton (ra-th biexciton) state fre- 



respectively, with V X j^ x (uj) = V X jf X (uj,uj). According to 
Eq. Q, the II rate depends on intensive parameters and 
therefore has finite value in the bulk limit. In contrast, 
the Auger rate has uncompensated volume prefactor and 
vanishes as V~ 2 in the bulk limit. The intraband cooling 
rates in the kinetic equations are calculated using Ohmic 
spectral density assuming no phonon bottleneck present 
in NCs. The total QE is evaluated through numerical so- 
lution of the kinetic equations with photogenerated pop- 
ulations used as the initial conditions. 

To calculate the electron and hole wavefunctions in 
spherically symmetric PbSe NCs and bulk, we use the 
effective mass formalism due to Kang and Wise [26] . 
The exciton and biexciton wavefunctions are introduced 
as uncorrelated configurations of the electron and hole 
wavefunctions to evaluate the transition dipoles and the 
interband Coulomb matrix elements (see Supplemental 
Material). The exciton (biexciton) dephasing rate is set 
to 50 meV (100 meV) which is order of magnitude consis- 
tent with Ref. [20 . To account for degeneracy of the exci- 
ton (biexciton) states associated with four equivalent L- 
valleys, the exciton (biexciton) population is multiplied 
by factor 4 (16). Introducing energy independent effec- 
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FIG. 1: Photogeneration QE vs pump energy for PbSe NCs 
of diameter d — 5,8,15,45 nm and bulk. Solid lines mark 
QE calculated using general expressions for populations from 
Ref. [24]. Dash marks approximate results (d=8, 15 nm NCs 
and bulk) due to Eq. The inset shows biexciton quantum 
yield associated with each photogeneration pathway for d = 
15 nm. 

tive electron-phonon coupling in the high DOS region, 
the intraband cooling rates are obtained by fitting the to- 
tal relaxation time, r p h, to reproduce the experimentally 
observed values in the range 0.5 < r p h < 5.0 ps [T3l [T9]. 
The Monte Carlo sampling is used to handle ensemble of 
NCs with 5% size distribution. 

Fig. [I] shows monotonic increase of the calculated pho- 
togeneration QE in the NC ensemble to the bulk values as 
the NC mean diameter increases. To explain this trent, 
we first look at the contributions of three photogenera- 
tion pathways (Eq. (|2|) to the biexciton quantum yield, 
Vxx = QE — 1. According to the inset to Fig. [T] the con- 
tribution of the last pathway (green line) is negligible. 
The other two become equal above the photogeneration 
energy threshold, ~ 1.5 eV. Taking into account reso- 
nant nature of the denominator in the first two terms 
of Eq. ([2]), the integral convolutions can be calculated. 
This results in the first and second term identical con- 
tributions to rj xx as observed in Fig. [IJ Hence, the total 
biexciton population acquires very simple dependence 

N xx (uj prn ) = kii(cj prn )N x (cj pm )/^ e ff, (6) 

on the II rate, /c//, photogenerated exciton population, 
N x , and the effective exciton-biexciton dephasing rate, 
Jeff- As a result, the photogeneration QE becomes de- 
pendent only on the II rate and effective dephasing rate, 
QE = 1 + kii/j e ff. To check this relation, we calcu- 
lated QE using II rate described by Eq.j4| (Fig. |5| with 
hjeff — 50 meV. According to Fig. fTjthe latter QE 
(dashed lines), well reproduces the trends of the "exact" 
calculations (solid lines). 

The analysis above leads to an important conclusion 
that in the absence of interference destroyed by high DOS 
and the interband dephasing processes, the photogenera- 
tion becomes resonant process. Furthermore, the biexci- 
ton generation described by Eq. (|6| can be interpreted as 



FIG. 2: Calculated (a) intensive biexciton DOS, (inset) effec- 
tive Coulomb interaction, and (b) II rate vs energy for PbSe 
NCs of diameter d — 8, 15 nm and bulk. 



a single II event of the exciton population occurring on 
the dephasing timescale, %ff Previous studies of photo- 
generation processes [5j [14] , implicitly assumed construc- 
tive interference within the second and the third path- 
ways. The assumption led to significant overestimation of 
photogeneration QE. Established linear dependence be- 
tween the photogenerated QE and II rate (Fig.[2jb)) is a 
key to explain the QE variation with NC size observed in 
Fig. [I] According to the inset to Fig. |2|a), the effective 
Coulomb interaction (Eq. (|3|) entering II rate (Eq. Q) 
is enhanced in NCs compared to the bulk limit just by a 
factor of two. In contrast, our calculations show that the 
biexciton DOS (Fig.[2|a)) reduction in NCs due to quan- 
tum confinement fully overplays the effective Coulomb 
enhancement. The dominant role of the DOS in II rate 
leads to the monotonic increase of QE in NC with their 
size increasing as observed in Fig. [T] This has also been 
found in Ref. [13] 

In contrast to II rate (Eq. p|) , the volume scaling of 
Auger recombination rate (Eq. (pi) makes this competing 
process negligible as the NC diameter increases. There- 
fore, the QE associated with the population relaxation 
is fully determined by the interplay between compara- 
ble II and cooling, r p h times. As a result, the total (i.e., 
photogeneration and population relaxation) QE has the 
same trend with the NCs size variation (Fig. [3| as II rate 
(Fig |2^b)) and never exceeds the bulk values. For com- 
parison, the experimental data for PbSe NCs [TT] and 
the bulk [13 are also shown in the plot. Solid lines 
in Fig. [3] show QE calculated for the intraband cool- 
ing time r p h = 1.0 ps which provides the best fit for 
the bulk. [13] Generally, the theory reproduces the ex- 
perimental trends in bulk (also supporting the adopted 
model for the phonon assisted cooling). Note that for 
tiWph ^ 3.6 eV where the tri-exciton generation appears, 
the application of our model becomes limited. To get 
better agreement with experiment for NC of d = 8 nm, 
we had to increase the cooling time up to r p h = 2 ps 
(blue dash). The increase of r p h can be attributed to 
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FIG. 3: Total QE as a function of pump energy calculated for 
PbSe NCs of diameter d = 8, 15 nm and bulk. Solid (dashed) 
lines present calculations with r p h = 1 ps (r p h = 2 ps). For 
comparison, diamonds and dots mark experimental data from 
Refs. [11] and [13], respectively. 



the quantum confinement induced increase in the level 
spacing [27] . The photogeneration contribution (Fig. [TJ 
to the total QE (Fig. [3| is small (~ 5%). This follows 
from the the fact that the photogeneration is a single II 
event occurring on short (< 100 fs) dephasing timescale. 
In contrast, during the population relaxation multiple II 
events take place on much longer (> 1 ps) timescale and 
provide major contribution to QE. 

To conclude, use of IESM parameterized by the ef- 
fective mass Hamiltonian allowed us to investigate the 
mechanisms involved in the photogeneration and popu- 
lation relaxation processes contributing to CM in PbSe 
NCs and bulk. We argue that direct photogeneration has 
the impact ionization nature and makes weak contribu- 
tion to the total QE. The use of intensive DOS and effec- 
tive Coulomb interaction further allowed for a rigorous 
analysis of the QE scaling with NCs size and comparison 
with the well denned bulk limit. We have found that QE 
in NCs plotted on absolute energy scale does not exceed 
that in bulk. This is in agreement with reported experi- 
mental data and confirms that the quantum-confmement- 
induced reduction in the biexciton DOS makes dominant 
contribution to QE. Finally note that the power per- 
formance of photovoltaic devices depends on the band 
gap scaling, and even for the case of weak Coulomb en- 
hancement, the NCs are potentially good candidates for 
such applications. Increase in the effective Coulomb in- 
teraction and biexciton DOS through the nanoscale het- 



erostructuring and surface functionalization should fur- 
ther improve their prospective performance. 
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